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' Abstract. We present full QCD and SUSY-QCD corrections to neutralino pair annihilation into 

\Q \ quark-antiquark pairs. We show results and numerically evaluate their impact on the neutralino 

annihilation cross section and the relic density in scenarios with non-universal gaugino masses. By 
comparing to the relic density obtained from a pure leading order calculation, we demonstrate that 
the corrections strongly influence the extraction of SUSY mass parameters from cosmological data. 
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INTRODUCTION 



The cosmological observations such as the mission of the Wilkinson Microwave 
Anisotropy Probe (WMAP) which indicate the existence of Cold Dark Matter (CDM) in 
the universe, provide first hints of physics beyond the Standard Model (SM) as all SM 
particles are excluded to be the dark matter candidate. On top of that the measurement of 
the relic density of dark matter has become so precise that it turned into a very powerful 
tool to test models beyond the Standard Model which contain a suitable candidate for 
dark matter. A combination of the five-year measurement of the cosmic microwave 
background by the WMAP mission with supernova and baryonic acoustic oscillation 
data yields the narrow 2a-interval for the relic density of dark matter JH] 



0.1097 <Q C DMh 2 < 0.1165, (1) 

where h denotes the present Hubble expansion rate Hq in units of 100 km s _1 Mpc -1 . 

The relic density of the dark matter particle % is directly proportional to its present 
number density hq. It is obtained by solving the Boltzmann equation describing the time 
evolution of the number density 

= - 3Hn% - (a ann v) {n\ - n 2 eq ) , (2) 

where the first term on the right-hand side corresponds to a dilution due to the expan- 
sion of the universe, and the second term corresponds to a decrease due to annihilations 
and co-annihilations of the relic particle into SM particles [2]. The Boltzmann equation 



TABLE 1. High-scale parameters together with the corresponding neutralino relic density, 
the contribution from top quark-antiquark final states to the annihilation cross section obtained 
with micrOMEGAs 2 . 1 [7] for our two selected non-universal gaugino mass scenarios. 
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can be integrated numerically and the energy density of the dark matter particle is ap- 
proximately inversely proportional to the thermally averaged annihilation cross-section 

(<7 ann v). 

Here, we use the precise experimental measurement of dark matter's relic density to 
constrain parameter space of the Minimal Super symmetric Standard Model (MSSM). In 
order to keep up with current and future experimental improvements, one has to reduce 
the theoretical uncertainty involved in the analysis, which we do by including the higher- 
order supersymmetric QCD corrections to annihilation processes with quark final states 
(for details see Ref. yfl). 

In many scenarios of the MSSM, the dark matter candidate is the lightest neutralino. 
In our analysis, we concentrate on scenarios that are derived from minimal supergravity 
(mSUGRA) models, where heavy quark final states dominate the annihilation. We 
choose to relax the gaugino mass unification which allows to investigate scenarios 
where the cross-section is dominated by the exchanges of the Z°-boson or the /-channel 
scalar quark exchanges as opposed to the Higgs-boson exchange in mSUGRA (see 
Refs. iSi). 

The soft SUSY-breaking Lagrangian containing gaugino mass terms for the 5-ino, 
V7-ino, and gluino has the following form 

^soft C -^(m 1 BB+M 2 WW +M 3 gg + h.c.). (3) 

Relaxing the universality of gaugino masses at the GUT scale, the values of M\, M 2 , 
and M3 at the unification scale can be considered as independent parameters. The 
sfermion and Higgs sectors are parametrized by the usual mSUGRA parameters mo, 
m\ji, Aq, tan/3, and sgn(/i) at the unification scale. Here, we adopt a commonly used 
parametrization and introduce the dimensionless parameters 

Xl = W 2 and X3 = W 2 ' (4) 

which will be used together with the W-ino mass parameter M 2 to describe the gaugino 
sector. The case x\ = x$ = 1 reproduces the mSUGRA model. 



RESULTS 

We investigate two scenarios with non-universal gaugino masses where the relevant 
parameters are given in Tab. [TJ The common feature of both scenarios is a lower gluino 
mass M3 which leads through the renormalization group evolution and electroweak 
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FIGURE 1. The contributions of the different diagrams to the annihilation cross section of a neutralino 
pair into top quark-antiquark pairs (left) and the effect of the radiative corrections on the annihilation 
cross section (right) as a function of the center-of-momentum energy p cm for the parameter points I (top) 
and II (bottom). The shaded area indicates the velocity distribution of the neutralino at the freeze-out 
temperature in arbitrary units. 



symmetry breaking to an increased higgsino fraction of the neutralino. An increased 
higgsino fraction enhances the Z°-boson exchange as in our scenario I (see top row 
of Fig. [T]) and also the squark exchange is larger due to enhanced Yukawa couplings. 
Moreover, the mass of the scalar tops can be lowered by decreasing the scalar mass 
parameter mo and using the fact that the gluino mass parameter is independent. Very 
light stops enhance their /-channel exchange and lead to dominant quark final states as 
in our scenario II (see bottom row of Fig.Q]). 

The effects of SUSY-QCD radiative corrections on the relic density (blue) is shown 
in Fig.[2]for both scenarios from Tab.Q]and compared to the tree-level result (green) and 
to the approximation implemented in the public computer code micrOMEGAs 2 . 1 
(red). For both scenarios there are important contributions from squark and Z° -boson 
exchanges due to the large destructive interference between the two channels. In general, 
the SUSY-QCD corrections are between 15%-20% and are larger than the experimental 
uncertainty. Therefore, distinct bands appear in scans around our parameter points I 
and II (see Fig. [2]). This demonstrates that the SUSY-QCD corrections are a necessary 





FIGURE 2. Cosmologically favored regions in the mo-M2-plane (top and bottom left) for the parameter 
points I and II. We show the regions that satisfy the constraint from Eq. ([TJ for our tree-level calculation 
(green), the calculation implemented in micrOMEGAs (red), and our calculation including the full SUSY- 
QCD corrections (blue). The prediction of the neutralino relic density iicDM^ 2 including the tree-level 
(green dash-dotted) cross section, the approximation included in micrOMEGAs (red dashed), and the full 
one-loop SUSY-QCD corrected cross-section (blue solid) as a function of either the high-scale gaugino 
mass parameter Mj_ for fixed mo = 300 GeV (top right) or the stop mass for a fixed neutralino mass 
my, = 235.6 GeV (bottom right). 



ingredient when identifying favored regions of parameter space or when extracting the 
underlying SUSY-breaking parameters from cosmological data. 
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